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Recently, charge radii and ground-state electromagnetic moments of Li and Be isotopes were 
measured precisely. We have performed large-scale ah initio no-core shell model calculations for 
these isotopes using high-precision nucleon-nucleon potentials. The isotopic trends of our computed 
charge radii, quadrupole and magnetic-dipole moments are in good agreement with experimental 
results with the exception of the ^^Li charge radius. The magnetic moments are in particular well 
described whereas the absolute magnitudes of the quadrupole moments are about 10% too small. 
The small magnitude of the ®Li quadrupole moment is reproduced, and with the CD-Bonn NN 
potential also its correct sign. 

PACS numbers: 21.60.De, 21.10.Ft, 21.10.Ky, 21.30.Fe, 27.20.+n 
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Recent developments of both experimental and theo- 
retical techniques have allowed for very precise measure- 
ments of charge radii and ground-state electromagnetic 
moments of exotic isotopes [3, H, In particular, 
charge radii can be determined through measurements 
of isotope shifts in certain atomic transitions using high- 
precision laser spectroscopy at isotope separation facili- 
ties [l|, Q. Theoretical calculations of electron correla- 
tions, as well as relativistic and QED corrections, can 
be performed [H, @] to yield very small uncertainties in 
the extracted results. Electric quadrupole and magnetic 
dipole moments can be determined using an experimental 
method that is based on the nuclear magnetic resonance 
technique 0, 

These observables reflect, in different ways, the evolv- 
ing nuclear structure along the isotopic chains. Fascinat- 
ing trends with varying N/Z ratio have been revealed [t^. 
It is a true challenge for theoretical methods to compute 
these observables, and to reproduce all of the observed 
trends simultaneously. In particular, for ab initio many- 
body methods this type of study allows to test properties 
of the high-precision nuclear Hamiltonians that are used 
as the single input to the calculations. But these calcu- 
lations also constitute a critical test of the limitations of 
the many-body method itself and of the operators that 
are used to compute the matrix elements. 

Efforts to compute some of these observables with mi- 
croscopic approaches have been performed using different 
realistic and semi-realistic interactions [H, 0, [iqI.Ii1| . Con- 
cerning the ab initio no-core shell model (NCSM) used in 
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this study, increased computational capabilities and im- 
proved algorithms have led to the opportunity to reach 
significantly extended model spaces. In addition, studies 
of realistic nuclear Hamiltonians have led some authors 
to explore the extent to which effects of multi-nucleon 
forces can be absorbed by non-local terms in the NN 
interaction [l2|, The use of these NN interactions 

allows, to some extent, to study three-body interaction 
effects while still maximizing the size of the model space. 

In concert, these developments can be regarded as 
a strong motivation to repeat and extend some earlier 
NCSM studies 0, [lH- In this Rapid Communication 
we calculate the charge radii and electromagnetic mo- 
ments of the A < 11 chains of Li and Be isotopes. We 
compare the performance of two very different NN inter- 
actions: (i) the CD-Bonn 2000 interaction (CDB2k) [3, 
that is a charge-dependent NN interaction based on one- 
boson exchange; and (2) the INOY IS-M [lH that is a 
phenomenological interaction for which non-locality was 
introduced in certain partial waves so that the binding 
energies of and ^He are described correctly [l^- All 
calculations are performed up to very large model spaces 
and efforts are made to quantify the rates of convergence 
of the observables. 

The NCSM method has been described in great detail 
in several papers, see e.g., Refs. [l8|, [l9|, [20|. Here, we 
just outline the approach as it is applied in the present 
study. We start from the intrinsic two-body Hamiltonian 
for the A-nucleon system Ha = T^ei + V, where T^ei is the 
relative kinetic energy and V is the sum of two-body nu- 
clear and Coulomb interactions. We solve the many-body 
problem in a large but finite harmonic-oscillator (HO) 
basis truncated by a chosen maximal total HO energy of 
the A-nucleon system. If we use realistic nuclear inter- 
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actions such as CDB2k or INOY, that generate strong 
short-range correlations, it is necessary to compute a 
model-space dependent effective Hamiltonian. For this 
purpose, we perform a unitary transformation [2l|, [22[ of 
the Hamiltonian, which accommodates the short-range 
correlations. In general, the transformed Hamiltonian 
is an A-body operator. Our simplest, yet nontrivial, 
approximation is to develop a two-particle cluster effec- 
tive Hamiltonian, while the next improvement is to in- 
clude three-particle clusters, and so on. The effective 
interaction is then obtained from the decoupling condi- 
tion between the model space and the excluded space for 
the two-nucleon transformed Hamiltonian. Details of the 
procedure are described in Refs. fiol. [23^ . 
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FIG. 1: (Color online) /iQ- dependence for ^Li ground-state 
observables calculated with the INOY interaction. Each curve 
corresponds to a particular model space represented by the 
parameter A^max (see text for details). 

In the present application we aim to describe nuclear 
observables that, in some cases, are sensitive to long- 
range properties of the wave function. It is therefore 
essential to include as many terms as possible in the ex- 
pansion of the total wave function. For this reason, we 
limit ourselves to the use of realistic two-nucleon (A/TV) 



interactions and we only present results obtained with 
the two-body cluster approximation. 

Our calculations are not variational, so higher order 
terms may contribute with either sign to total binding. 
Hence, evaluating the dependence on the basis-space pa- 
rameters help calibrate our convergence. Once the ef- 
fective interaction is derived, we diagonalize the effec- 
tive Hamiltonian in a Slater determinant HO basis that 
spans a complete A/'max^^ space. This is a highly non- 
trivial problem because of the very large dimensions we 
encounter. To solve this problem, we have used a special- 
ized version of the shell model code Antoine Q , adapted 
to the NCSM [25[. The first step in our study is a search 
for the optimal HO frequencies. The converged results 
should not depend on the HO frequency. In practice 
however, due to the cluster approximation of the effective 
interaction, our results will be sensitive to the choice of 
hCl. Furthermore, by construction the effective interac- 
tions depend on the size of the model space, A/'max, and 
on the number of nucleons, A. In order to investigate 
these dependences we have performed calculations for a 
large series of frequencies. As an example, the results 
from this study for the ^Li ground-state binding energy 
and electromagnetic moments are presented in Fig. [1] for 
the INOY interaction. 

We are looking for the region in which the dependence 
on HQ is the smallest; and we select this frequency (from 
the calculation in the largest model space) as a starting 
point for our detailed investigation of ground-state ob- 
servables. The selected optimal HO frequencies, HQ a, 
are presented in Table U 



TABLE I: Selected optimal HO frequencies, HQ a, (in [MeV]). 
Also shown is the maximum model space (represented by 



A^max) that 


was reached for each mass 


number A. 
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^The choice fiOA=12 MeV is used for ^^Li. 

In addition, our results will depend on the size of the 
model space represented by the parameter A^max- Our 
approach is to fix the HO frequency, at the value HQa 
specified in Table [H and observe the convergence with 
increasing model space. At the same time, we investigate 
the A/'max-dependence at HQa ± 1 MeV. 

Using these studies we can gauge the degree of conver- 
gence of our final results. For a given observable O the 
total error AO is estimated by 



AO 



AO^ 



AO^, 



(1) 



where AOat and AOq represent the observed rates of 
convergence with respect to model-space size and HO 
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frequency, respectively. We choose AOn simply as the 
difference between the numerical results in the largest 
and next-to-largest model spaces using the optimal HO 
frequency. To determine AOq we consider the triplet 
of results computed at HQa and HQa i 1 MeV in the 
largest model space for which all three are available. We 
then define the quantity AOq as the difference between 
the maximum and minimum value of this triplet of re- 
sults. The total error AO is presented in the tables of 
this paper. Ground-state binding energies {E) for Li iso- 
topes are summarized in Table [III We note that for the 

TABLE IL Ground- state binding energies (E) for Li iso- 
topes. The INOY results are extrapolated energies (with 
convergence- rate estimates), while the actual computed re- 
sults from the largest model space are given in square brack- 
ets. The experimental results are from i27|, i28i] . 





GDB2k 


E [MeV] 
INOY 




Exp 


^Li 


29.07(41) 


32.33(19) 


[32.07] 


31.99 


^Li 


35.56(23) 


39.62(40) 


[38.89] 


39.24 


«Li 


35.82(22) 


41.27(51) 


[39.94] 


41.28 


^Li 


37.88(82) 


45.86(74) 


[42.30] 


45.34 


i^Li 


37.72(45) 


42.50(95)" 


[40.44] 


45.72(1) 



*The exponential convergence rate is not fully reached. 



INOY NN potential, the ground-state energy convergence 
is very uniform with respect to the HO frequency with 
systematic changes with A^max • The convergence with in- 
creasing A^max is quite evident in particular for ^4 = 6 — 9, 
see e.g. the top panel of Fig. [U and we are able to ex- 
trapolate assuming an exponential dependence on A'max 
as £^(iVmax) = Eoo + a exp(-67Vmax). The same type 
of extrapolation was successfully performed for the He 
isotopes [29]. Therefore, extrapolated INOY energies are 
given in Table [TT] together with computed results from 
the largest model space (in square brackets). 

It also deserves to be mentioned that we can, in prin- 
ciple, have different fir^-dependence for different observ- 
ables. In particular since we're using effective Hamiltoni- 
ans to compute the eigenenergies and eigenfunctions but 
are employing bare operators to evaluate other observ- 
ables. In practice, however, the choice of HQ a is relatively 
stable. The exception to this general rule is observed for 
those observables that depend critically on the size of 
the system. In particular, for large systems we some- 
times observe a different fil^-dependence for the electric 
quadrupole moment and for the radius. In this case the 
convergence-rate estimate should be viewed as provid- 
ing the magnitude of variations around the maximum 
model space reached. The effective-operator formalism 
for general one- and two-body operators was developed 
and studied in Ref. [sO] , finding weak renormalization of 
long-range operators. 

To complete the present survey we have collected re- 
sults for odd-A Be isotopes from earlier papers. The 



''Be results are partly from Ref. [ij while ^'^^Be results 
are from Ref. [l5|. However, additional calculations were 
performed in the present study to provide convergence 
rate estimates. In addition, new calculations have been 
performed for the ^^Be ground state. 

Ground-state quadrupole and magnetic dipole mo- 
ments have been calculated in the impulse approxima- 
tion using bare operators. Two-body meson-exchange 
currents have been shown to increase isovector magnetic 
moments [3l|, |32| but are not included here. The A^max- 
dependence of Q and /i is presented in Fig. [2] for odd 
Li isotopes using the CDB2k interaction. We note that 
the magnetic moment is usually very well converged in 
the NCSM. The quadrupole moment is sensitive to the 
length scale of the NCSM basis and for smaller model 
spaces it tends to increase with decreasing HO frequency. 
Final results for all isotopes are presented in Table Hill to- 

TABLE III: Ground-state quadrupole moments (Q), mag- 
netic dipole moments (/x), and charge radii (rc) for Li and 
Be isotopes. Theoretical results for ^Be and ^'^^Be are 
partly from Refs. p^. [TH . The experimental results are from 
Refs. S 05 [23, [23 for electromagnetic moments, and from 
Refs. [il, 13, @] for radii. 







Q [e fm^] 






CDB2k 


INOY 


Exp 


^Li 


-0.066(40) 


+0.080(19) 


-0.0806(6) 


^Li 


-3.20(22) 


-2.79(17) 


-4.00(3) 


«Li 


+2.78(12) 


+2.55(12) 


+3.14(2) 


^Li 


-2.66(22) 


-2.30(12) 


-3.06(2) 




-2.81(27) 


-2.32(13) 


-3.33(5) 




-5.50(48) 


-4.68(28) 






+4.12(26) 


+3.67(23) 


+5.288(38) 






/i [/iTv] 






CDB2k 


INOY 


Exp 


^Li 


+0.843(5) 


+0.843(2) 


+0.822 


^Li 


+3.01(2) 


+3.02(2) 


+3.256 


«Li 


+ 1.24(6) 


+ 1.42(4) 


+ 1.654 


^Li 


+2.89(2) 


+2.98(5) 


+3.437 




+3.56(4) 


+3.54(4) 


+3.671(1) 


^Be 


-1.14(1) 


-1.15(1) 


-1.3995(5) 


^Be 


-1.22(9) 


-1.06(6) 


-1.1774 


i^Be 


-1.55(6) 


-1.47(3) 


-1.6813(5) 






rc [fm] 






CDB2k 


INOY 


Exp 


^Li 


2.40(6) 


2.29(4) 


2.540(28) 


^Li 


2.36(7) 


2.20(5) 


2.390(30) 


«Li 


2.31(8) 


2.16(5) 


2.281(32) 


^Li 


2.25(10) 


2.07(6) 


2.185(33) 




2.26(13) 


2.06(6) 


2.426(34) 


^Be 


2.56(10) 


2.39(5) 


2.647(17) 


^Be 


2.41(11) 


2.22(6) 


2.519(12) 


lOBe 


2.34(9) 


2.16(6) 


2.357(18) 


i^Be 


2.37(11) 


2.19(7) 


2.463(16) 



gether with recent and very precise experimental results 
using the nuclear magnetic resonance technique 0, 0]. 
We note that the two different interactions used in this 
study give very similar isotopic trends but with a consis- 
tently smaller magnitude for the INOY interaction. This 
observation is connected to the anomalously large nuclear 
density generated by this interaction found already in 
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^He calculations [i^,!!^. For ^Li the quadrupole moment 
is basically converged and the magnitude agrees well with 
the experimental value. However, for the larger systems 
the calculated magnitudes are consistently > 10% too 
small as compared to the experimental results. This re- 
sult is expected from working in a finite HO basis and we 
observe that the magnitude of the quadrupole moment is 
steadily increasing with increasing model space, as shown 
in Fig. [2 

The point-nucleon radii {vpt-p and Vpt-n) can be com- 
puted from the NCSM wave functions with complete re- 
moval of spurious CM motion. We have translated these 
results into a rms charge radius {vc) by foldin g in pro- 
ton {R^) and neutron {R^) rms charge radii [SJ and 
adding the Darwin-Foldy term [35[ according to the re- 
lation {rl) = {rl,_^) + (i?^) + N{Rl)/Z + my{4M^c^). 
These results are compared to recent high-precision laser 
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FIG. 2: (Color online) A^m a cc- dependence of ground-state elec- 
tromagnetic moments and charge radius for odd Li isotopes 
calculated with the CDB2k interaction. Each curve corre- 
sponds to a fixed HO frequency. In particular, thick (thin) 
lines correspond to HQ a (h^A i 1 MeV) for each isotope. 

spectroscopy measurements di H in Table lUll The ex- 
perimental results are obtained using previously mea- 



sured ^Li and ^Be radii as references. The convergence 
of the calculated radius with increasing A/'max is shown in 
Fig. [2] for odd Li isotopes using the CDB2k interaction. 
We also conclude that the larger binding energy obtained 
with the INOY interaction can partially be explained by 
an abnormally large nuclear density and consequently a 
charge radius that is too small. For CDB2k, the radius re- 
sults are rather well converged and show good agreement 
with experiment. However, the known halo structure of 
the ^^Li ground-state is not reproduced in the limited 
HO model space of the NCSM. 

In summary we have computed charge radii and elec- 
tomagnetic moments of Li and Be isotopes using two 
different, high-precision nuclear Hamiltonians within the 
ab initio NCSM. In Fig. [3] we compare the calculated 
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FIG. 3: (Color online) NCSM calculated electric quadrupole 
moments, magnetic dipole moments, and charge radii of Li 
isotopes compared with experimental results. See also Ta- 
ble ini 

and experimental trends for these observables for the Li 
chain of isotopes. With the exception of the radius of the 
^^Li halo ground-state we find a very good agreement be- 
tween NCSM results and recent experiments. The overall 
trends of all observables are well reproduced. Magnetic 
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dipole moments are usually characterized by very good 
convergence properties in the NCSM and we find a good 
agreement with the experimental values. Another suc- 
cess is the tiny quadrupole moment of ^Li that is known 
to pose a difficult task for most theoretical calculations. 
In particular, the general failure of three-body models 
for this observable has been blamed on missing anti- 
symmetrization of the valence nucleons and the nucleons 
in the alpha-core [s^. The NCSM correctly reproduces 
the very small value, but with CDB2k and INOY giving 
different signs. Simultaneously, the trend for the much 
larger moments of A = 7 — 11 is nicely reproduced. We 
note that the ratio Q(^^Li)/Q{^Li) is found to be very 
close to unity, as confirmed recently by very precise ex- 
perimental data [3]. This finding is obtained without a 
very accurate description of the dilute halo structure of 
-^^Li; a structural feature that we find would require an 
extension of the HO basis used in the standard NCSM. 
Still, the decrease of the charge radius of A = 6 — 9 iso- 



topes is reproduced, although the INOY interaction gives 
too high nuclear densities. 

We conclude by observing that the recent achievement 
of performing very precise measurements of ground-state 
properties of exotic isotopes proves to be a very valu- 
able tool in the quest for understanding the nuclear in- 
teraction, the forms of relevant operators, as well as the 
evolving structure of the nuclear many-body system. 
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